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Abstract The measurements of pulsar frequency sec- 
ond derivatives have shown that they are 10 2 — 10 6 times 
larger than expected for standard pulsar spin-down law, 
and are even negative for about half of pulsars. We ex- 
plain these paradoxical results on the basis of the statis- 
tical analysis of the rotational parameters is, v and v of 
the subset of 295 pulsars taken mostly from the ATNF 
database. We have found a strong correlation between v 
and v for both i) > and v < , as well as between v and 
v . We interpret these dependencies as evolutionary ones 
due to v being nearly proportional to the pulsars' age. 
The derived statistical relations as well as "anomalous" 
values of v are well described by assuming the long-time 
variations of the spin-down rate. The pulsar frequency 
evolution, therefore, consists of secular change of v ev (i), 
v ev {t) and u ev (t) according to the power law with n w 5, 
the irregularities, observed within a timespan as a tim- 
ing noise, and the variations on the timescale larger than 
that timespan - several tens of years. 
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1 Introduction 

The spin-down of radio pulsars is due to the conver- 
sion of their rotation energy into emission. According 
to the "classical" approach, their rotational frequencies 
v evolve obeying the spin-down law v = —Ki/ n , where 
K is a positive constant that depends on the magnetic 
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dipole moment and the moment of inertia of the neutron 
star, and n is the braking index. The latter can be de- 
termined observationally from measurements of u, v and 
v as n = vv/v 2 . For a simple vacuum dipole model of 
pulsar magnetosphere n — 3; the pulsar wind decreases 
this value to n = 1; for mu ltipole magnetic field n > 5 
(|Manchester fc Tavlo dfl977h . At the same time the mea- 
surements of pulsar frequency second derivatives v have 
shown that their values are much larger than expected 
for standard spin-down law and are even negative for 
about half of all pulsars. The corresponding braking in- 
dices range from -10 6 to 10 6 dD'Alessandro et ail 1 19931: 
IChukwuddll993l : iHobbs et al.ll2004D . " 

It was found that the significant correlations between 
\i>\ (\P\) and v (P) demonstrate the increase of the abso- 
lute values of the v and P second derivatives for younger 
(with faster slow-down) pulsars dCordes &, Downslll985t 
lArzoumanian et al.lll994j : lLvndfl999h . The anomalously 
high and negative values of v and n may be interpreted 
as a result of low-frequency terms of the "timing noise" 
- a complex variations of pulsars rotat ional phase within 
a timespan ( D'Alessandro et aT1ll993l ). 

It is clear that the timespan of observations is by no 
means intrinsic to the pulsar physics. Indeed, the varia- 
tions of rotational parameters may take place on larger 
timescales as well. However, the timescale of observations 
naturally divides it into two separate classes of manifes- 
tations - the residuals in respect to the best fit for the 
timing solution (the "timing noise") and the systematic 
shift of the best fit coefficients (i.e. in the measured val- 
ues of u, i)) relative to some mean or expected value 
from the model. The latter effect may be called the "large 
timescale timing noise" . 

Up to date we know nearly 200 pulsars for which the 
timespan of observations is greater than 20 years, and 
the v alues of their v st ill turn out to be anomalously 
large (jHobbs et al.ll2004D . 

For example, for the PSR B 1706- 16 pulsar, variations 
of v with an amplitude of 10 -24 s -3 have b een detected 



on a several years timescale (see Fig. 7 in ( Hobbs et al.l 
l2004f )). with the value of v depending on the time interval 
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168 pulsars 
r = 0.90; SL : 
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127 pulsars 
r = 0.85;SL 



J (-v) 
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nearly linear laws, however they seem to be not exactly 
symmetric relative to v = Q. We divided both branches 
into 6 intervals of i>, computed mean values and its stan- 
dard deviations of v± in each and rejected the hypothe- 
sis of their symmetry with a 0.04 significance level. The 
absolute values of i> + are systematically larger than cor- 
responding V- (the difference is positive in 5 intervals of 
6). Also, the difference of analytical fits to branches is 



positive over —10 



-10 



-15 



interval of v. These 



Fig. 1 The v — v diagram for 295 pulsars. T he figure shows 
the pulsars from the work (Hobbs et al. 2004) as circles, and 
the objects measured by other groups as squares. Open sym- 
bols represent the pulsars associated with supernova rem- 
nants, and therefore - relatively young ones. Analytical fits 
for both positive and negative branches are shown as solid 
lines. Measurment errors shown as an error bars. 



selected. However, the fit over the entire 25 year times- 
pan gives a value of v = 3.8 TO -25 s -3 with a few percent 



accuracy (which leads to a braking index « 2.7 ■ 10 3 ). 

In the current work we provide observational evi 
dence of the nonmonotonic evolution of pulsars on timescales^^^ any d " e p endenc " e " on 7o7r c 7 reflectlthe'de- 
larger than the typical contemporary timespan of obser- pendence on pulsar age 



are the arguments in favour of a small positive asymme- 
try of branches. 

We found an obvious correlation of v with the char- 
acteristic age T ch = -\K (r = 0.96, v ~ r^ 11 ^ 02 ). 
These parameters are nearly proportional, which leads 
to a significant correlation of r c h both with i) (r = 0.85 
for the positive branch and r — 0.75 for the negative one) 
and with n (r = 0.75 and r = 0.76 correspondingly). 

The correlation s found are fully consis t ent with the 

result s publis h ed in Cordes fc Downsj (1 19851): Arzoumanian et all 
(|1994D ; iLvnel (|l999f ). as well as in lUrama et all (j2006t ). 
However, the branches with v > and v < in those 
works were not analyzed separately from each other (not 
as 

Young pulsars confidently associated with supernova 
remnants are systematically shifted to the left in Fig. [T] 
(open symbols). The order of their physical ages roughly 
corresponds to that of their characteristic ages. This 



vations (tens of years), using the statistical analysis of 
measured v, v and v. We estimate the main parame- 
ters of such long timescale variations and discuss their 
possible relation to the low-frequency terms of the tim- 
ing noise. We have also derived the parameters of pulsar 
secular spin-down. 



2 Statistical analysis of the ensemble of pulsars 



Our statistical analysis is based on the assumption that 
numerous measurements of the pulsar frequency second 
derivatives reflect their evolution on the timescale larger 
than the duration of observations, and uses the parame- 
ters of 295 pulsars. 

From 389 objects of ATNF catalogue ( Manchester et all increasing of v implies a positive evolutionary value of v 

2005) with known v we have compiled a list of " ordinary" The characteristic timescale T of such variations must 
radio pulsars with P > 20 ms and P > 10~ 17 s s~ l , be much shorter than the pulsar life time and at the 
excluding recycled, anomalous and binary pulsars, and same time much larger than the timescale of the obser- 



The v — v diagram (Fig. [T]) may be interpreted as 
an evolutionary one. In other words, each pulsar during 
its evolution moves along the branches of this diagram 
while increasing the value of its v (which corresponds to 
the increase of its characteristic age). However, there is 
an obvious contradiction: for the negative branch, is, be- 
ing negative, may only decrease with time (since v is for- 
mally the derivative of u), and the motion along the neg- 
ative branch may only be backward! This contradiction 
is easily solved by assuming non-monotonic behaviour of 
v(t), which has an irregular component (Sv) along with 
the monotonic one (v e v), where the subscript "ev" marks 
the evolutionary value. In such interpretation the value 
of i) ev must be positive, which will lead to a positive 
asymmetry of the branches. Moreover, the evolutionary 



with relative accuracy of second derivative measurements 
better than 75%. I t has been appended with 26 pulsars 
from other sources (jD'Alessandro et aTlll993t IChukwudel 
1993). The parameters of all pulsars have been plotted 
on the v — v diagram (Fig. [l| . 

The basic result of the statistical analysis of this data 
is a significant correlation of v and v : both for 168 ob- 
jects with v > (correlation coefficient r w 0.90) and for 
127 objects with v < (r w 0.85). Both groups follow 



vations. As it evolves, a pulsar repeatedly changes sign 
of i), in a spiral-like motion from branch to branch, and 
spends roughly half its lifetime on each one. The asym- 
metry of the branches reflects the positive sign of v ev (t), 
and therefore, secular increase of v e v{t) (i-C- ah pulsars 
in their secular evolution move to the right on the v — v 
diagram). Systematic decrease of branches separation re- 
flects the decrease of the variations amplitude and/or the 
increase of its characteristic timescale. 
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Any well known non-monotonic variations of v(t), like 
glitches, microglitches, timing noise or precession, will 
manifest themselves in a similar way on t he v— i> diagram 
and l e ad to extremely h igh values of v ( Shemar fc Lvnd 
Il996l IStairs et al . 2000). However, their characteristic 
timcscales vary from weeks to years, and they are de- 
tected immediately. But here the variations on much 
larger timescales are discussed, and their study is possi- 
ble only statistically, assuming the ergodic behaviour of 
the ensemble of pulsars. 



Correlation coefficient 1-O.68, SL = 10" 



• pulsars with v > 
pulsars with V < 




-15.27 + 1.35 5.13 +0.34 



3 Non-monotonic variations of pulsar spin-down 
rate on large timescales 



Fig. 2 The v — v diagram for the pulsars with the measured 
second derivative. The filled symbols are objects with posi- 
tive v, the open ones - with negative. The behavior of both 
subsets is the same. The solid line represents the best fit, cor- 
responding to the n~5 braking index, the dotted ones - 1-cr 
complicated - periodic, quasi-periodic, or completely stochag an g e . 

tic. Generally, it may be described as a superposition 



Variations of the pulsar rotational frequency may be 



v{t) = v ev (t)+Sv(t), 



(1) 



where v ev (t) describes the secular evolution of pulsar pa- 
rameters and 5v(t) corresponds to irregular variations. 
Similar expressions describe the evolution of v and v af- 
ter a differentiation. The Si>(t) satisfies the obvious con- 
dition of zero mean value < Si)(t) >t~ over the times- 
pans larger than characteristic timescale of the varia- 
tions. The amplitude of the observed variations of v is 
re lated to th e dispersion of this process as asa = A$ = 

V< Wf >■ 

The second derivative values on the upper i) + and 
lower i>_ branches in Fig. [1] may be described as v±(t) = 
u ev {t) ± Aj>(t) for each pulsar. This equation describes 
some "average" pulsar, while the spread of points inside 
the branches reflects the variations of individual param- 
eters over the pulsar ensemble and reaches 4 orders of 
magnitude. 

The second derivative v is the only parameter signif- 
icantly influenced by the timing variations (see Section 
0]). Thus one can assume that the measured values of v 
and v may be considered to be evolutionary ones, v ev 
and v ev (since bv and Sv are small). 

Using relations described above, the secular behav- 
ior v(t) (or, v(v)) may be found by plotting the studied 
pulsar group onto the v — v diagram (Fig. [2]). The ob- 
jects with v > and v < are marked as filled and 
open circles, correspondingly. It is easily seen that the 
behavior of these two sub-groups is the same, which is 
in agreement with the smallness of the pulsar frequency 
variations in respect to the intrinsic scatter of v{v). How- 
ever, a strong correlation between v and v (r « 0.7) is 
seen, and 



-Cv n , 



(2) 



where C = 10- 15 26±1 - 38 and n = 5.13 ± 0.34. So, the 
secular evolution of the "average" pulsar is according to 
the "standard" spin-down law with n ~ 5! This result is 
very interesting on its own, especially since the v and v 



are always measured as independent values. The brak- 
ing index rj 5 may suggest the importance of multipole 
components of pulsar magnetic field, or the deviation of 
the a ngle between pulsar rotation al and dipolc axes from 
tt/2 ([Manchester fc Tavlorl [l977h . 

Note here, that the width of the fit on Fig. [5] is 
quite small - only about 1.5 orders of magnitude. If the 
spin-down is even approximately close to being described 
by the vacuum dipole model, then C should scale as 
{Bo sin a) 2 , where B is the polar field and a is the mag- 
netic inclination angle. But the range of (Bq sin a) 2 over 
the pulsar population is expected to be of many orders 
of magnitude. Therefore, Fig. [2] once again shows that 
simple vacuum dipole model is not adequate to observa- 
tions. The braking index n rj 5 allows the interpretation 
of C in terms of the quadrupole spin-down. It could be 
explained also in the frame of electric curre nt mechanism 
of pulsar spin-down ( Gurevich et al.lfl993h . 

As has been stated above, the dependency between 
i> and r c /j is consistent with a power law with a slope 
of —1.16 ±0.02. This value is significantly different from 
— 1.0 which would be the case when the v — v correlation 
is absent and roughly consists with a slope of 5.13 in 
the v{y) dependency (the values are different on a 2.5cr 
level). But n is strongly dependent on the i> — r slope 
value, so —1.16 seems to be more or less consistent with 
the measured value of n. 

From ([2]) we may easily determine the relation be- 
tween v ev and v as v ev — nC~ (— v) 2 ^ ~ , which is shown 
in Fig. [3] as a thick dashed line. The same relation may 
be also estimated directly by using the asymmetry of the 
branches seen on Fig. Q] as v(v) — \ {v+ + t'-), where i>± 
are defined in Fig. [T] Such estimation, while being very 
noisy, is positive in the — 10~ n -: — 10~ 15 s -2 range and 
agrees quantitatively with the previous one. 

The amplitude of the v oscillations, A- v , may be easily 
computed in a similar way, by using i) + and V- , as Aj> — 
\{v+ - V-). 
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:10 -2.11±0.15 ( . 1 



1 1 + 0.04 
80 + 0.01- 




Fig. 3 The D — v diagram with the simple variations model. 
The solid line is the amplitude A v of the frequency second 
derivative variations, the dashed - the secular term i> e „, and 
the dot-dashed lines are the envelopes of the oscillations u ev ± 
Ad with 1-cr ranges (dotted lines). Each pulsar spends the 
majority of its life time at or very near the envelopes. 



The behaviour of pulsars according to the derived re- 
lations is shown in Fig. [3l This simple variations model 
describes the observed branches, both positive and nega- 
tive, rather well. The absence of negative branch objects 



with 



v > 



-10- 



i.e. with T c h > 10 , we interpret 



as a prevalence of the second derivative's secular com- 
ponent over the varying one {Ay < i> ev ) in this region. 
Older pulsars begin to change the sign of v due to spin 
rate variations. 



4 Discussion and conclusions 

In general, it is impossible to estimate the amplitudes 
of the frequency and its first derivative variations A v 
and A v from the amplitude of the second derivative only 
(the knowledge of its complete power density spectrum 
is needed). However, if the spectral density is relatively 
localized and some characteristic timescale T of the vari- 
ations exists, it is possible to set some limits on it. A 
rough estimation is A v ~ A V T, A v ~ A V T and A v ~ 
ApT 2 . On a large timescale the variations can not lead 
to pulsar spin-up, so the variations of frequency first 
derivatives are much smaller than the secular ones, and 
Ac ~ A U T « i), so T « u/Ajj. So, for the -10~ 12 < 



v < 



-10 



-15 



range and corresponding values of A v 



T up ~ 10 11 s. Also, this characteristic timescale is obvi- 
ously larger than the timespan of observations, so 50 < 
T < 3 • 10 3 years. Assuming the constancy of T during 
the pulsar evolution and therefore the change of A v with 
time, we get A v ~ 10~ 3 -j- 10 -7 Hz. For such a model 
the pulsar frequency varies with the characteristic time 
of several hundred years and the amplitude from 10~ 3 
Hz for young objects to 10~ 7 Hz for older ones. 



The physical reasons of the discussed non-monotonic 
variations of the pulsar spin-down rate may be similar 
to the ones of the timing noise on a short timescale. Sev- 
eral processes had been pro posed for their explanation 
(|Cordes fc Greensteinlll98lT ) - from the collective effects 
in the neutron star superfluid core to the electric cur- 
rent fluctuations in the pulsar magnetosphere. Whether 
these processes are able to produce long timescale vari- 
ations is yet to be analyzed. On a short timescale, the 
pulsars show different timing behaviour. But on the long 
timescale their behavior seems to be alike. 

The argument in favor of the similarity between the 
discussed variations and the timing noise is the coinci- 
dence of the timing noise v amplitu de extrapolated ac - 
cording to its power spectrum slope ( Bavkal et al.lll999h 



to the time scale of hundreds of years, with the A u de- 
rived from our analysis for the same v 1 i.e. the same ages 
(see Fig. 3). 

At the same time, there are several low noise pulsars 
wit h large or negative y. Fo r 19 of 45 pulsars studied 
in ( D'Alessandro et al.lll995| ) the timing noise is nearly 



absent (RMS < 1 • 10 



consistent with the 



(Hobbs et alJ 12004 


), which are well 


v\ - v correlation ( 


Cordes & Downs! 



i>. This shows a possible difference between timing noise 
and the long timescale variations described above. 

In any case, the principal point is that all the pulsars 
evolve with long-term variations, and the timescale of 
such variations significantly exceeds several tens of years. 
That explains the anomalous values of the observed v 
and braking indices and gives reasonable values of the 
underlying secular spin-down parameters. 
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